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ABSTRACT 

We present results from three XMM-Newton observations of the M31 low mass X-ray 
binary XMMU J004314.4+410726.3 (Bo 158), spaced over 3 days in 2004, July. Bo 
158 was the first dipping LMXB to be discovered in M31. Periodic intensity dips were 
previously seen to occur on a 2.78-hr period, due to absorption in material that is 
raised out of the plane of the accretion disc. The report of these observations stated 
that the dip depth was anti-correlated with source intensity. In light of the 2004 XMM- 
Newton observations of Bo 158, we suggest that the dip variation is due to precession of 
the accretion disc. This is to be expected in LMXBs with a mass ratio < 0.3 (period 
< 4 hr), as the disc reaches the 3:1 resonance with the binary companion, causing 
elongation and precession of the disc. A smoothed particle hydrodynamics simulation 
of the disc in this system shows retrograde rotation of a disc warp on a period of ~11 
Porb, and prograde disc precession on a period of 29±1 P 0 rb- This is consistent with 
the observed variation in the depth of the dips. We find that the dipping behaviour is 
most likely to be modified by the disc precession, hence we predict that the dipping 
behaviour repeats on a 81±3 hr cycle. 

Key words: X-rays: general - Galaxies: individual: M31 - X-rays: binaries - Accre¬ 
tion: accretion discs - Methods: numerical 


1 INTRODUCTION 

Bo 158 is source number 158 in the catalo gue o f globu¬ 
lar clu sters that were identified in M31 bv iBattistini et all 
Jl987tl . Its X- ray counterpart was discovere d by the Einstein 
Observatory IITrinchieri fc FabbianaHl99li source number 
81), and is located at a = 00 h 43 m 14.2 s , 5 = 41° 07 / 26.3" 
foi Stefano et a~ • 2002 ). iTrudolvubov et al.l (: 2002j) identi¬ 
fied the X-ray source as a likely low mass X-ray binary 
(LMXB) with a neutron star primary; following their work, 
we will use the designation “Bo 158” to describe the X-ray 
sou rce h e re. 

iTrudolvubov et~al . 1200 2) report ~83% modulation in 
the 0.3-10 keV flux of Bo 158 on a 2.78 hour period during 
the ^60 ks 2002 January XMM-Newton observation. The 
modulation resembles the intensity dips seen in high incli¬ 
nation LMXBs due to photo-electric absorption of X-rays by 
materia l th at is raised above the body of the accretion disc 
llWhite fc Swankll982i) . The authors comment that the dip¬ 
ping is energy-independent, and discuss two possible mech¬ 
anisms: obscuration of the central X-ray source by highly 
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ionised material that scatters X-rays out of the line of sight, 
and the partial covering of an extended source by an opaque 
absorber. They also report ~30% dips in the 2000, June 
XMM-Newton lightcurve and ~50% dips in the 0.2-2.0 keV 
lightcurve of the 1991, June 26 ROSAT/PSPC observation. 
However, no significant dips were found in the 0.3-10 keV 
lightcurve of the 2001 June XMM-Newton observation; the 
authors placed a 2cr upper limit of 10% on the modulation. 

ITrudolvubov et alJ (120021) obtained fluxes in the 0.3-10 
keV band for the three XMM-Newton observations by fitting 
absorbed Comptonisation models to the combined spectra 
of the EPIC-pn, MOS1 and MOS2, and concluded that the 
depth of the intensity modulation was anti-correlated with 
the source luminosity. We present further XMM-Newton ob¬ 
servations and modelling results which suggest that the vari¬ 
ation in dipping behaviour may instead be due to precession 
in the accretion disc. Such behaviour is associated with the 
“superhump” phenomenon that is observed in interacting 
binaries where the mass rat io of the secondary to th e pri¬ 
mary is smaller than ~0.3 dWhitehurst fc Kinelll99ll) . Su¬ 
perhumps are briefly reviewed in Sect. |2] followed by de¬ 
tails of the observations and data analysis in Sect. [3 and 
our results in Sect. 21 Numerical modelling of the system 
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is discussed in Sect. 5; the system was simulated by a 3D 
smoothed particle hydrodynamics (SPH) code. We present 
our discussion in Sect. El and finally our conclusion in Sect. 7. 


Table 1. Journal of XMM-Newton observations of the M31 core. 
Al—A3 are available in the public archive, while P4—P6 are propri¬ 
etary observations. Two observations were made during 2004 July 
19; the first observation (a) started at 01:42:12, and the second 
(b) started at 13:11:22. 


2 SUPERHUMPS 


Superhumps were first identified in the superoutbursts of 
the SU UMa sub-class of cataclysmic variables. They are 
manifested as a periodic increase in the optical brightness 
on a p eriod that is a few perce nt longer than the orbital 
period ll VoalUl OtS IWarrierll l97.Tl) . SU UMas are a subclass 
of dwarf novae with short orbital periods (<2 hr) that ex¬ 
hibit particularly long, bright superoutbursts, separated by 
several outbursts that are typical of all dwarf novae; the su¬ 
peroutburst inte rvals are >5 times longer than the normal 
outbursts. irVogtlll980l ). ___ 

In the model proposed bv lOsakl jlQSfifl . a superoutburst 
is occurs when a normal outburst is enhanced by a tidal in¬ 
stability; this occurs when the outer disc reaches a 3:1 res¬ 
onance with the secondary. The disc is small at the start of 
the superoutburst cycle, well within the radius of tidal in¬ 
stability, and little angular momentum is removed by tidal 
interaction during the first normal outburst. Hence, the disc 
cannot accrete all of its mass onto the neutron star, and the 
size of the disc increases with su ccessi ve outbursts, until the 
3:1 resonance is reached llOsakillll989il . The additional tidal 
forces exerted on the disc by the secondary at this stage 
cause the disc to elongate and precess, and also greatly en¬ 
hance the loss of angular momentum, so that the disc con¬ 
tracts, and most of the disc material is snow-ploughed ont o 
the neutron star, causing the superoutburst dOsakil 1198 9*1. 
The disc precession is prograde in the rest frame, and the 
secondary repeats its motion with respect to the disc on 
the beat period between the orbital and precession periods, 
slightly longer than the orbital period. The secondary modu¬ 
lates the disc’s viscous dissipation on this period, giving rise 
to maxima in the optical lightcurve, known as superhumps. 

The requirement for the 3:1 resonance to fall within the 
disc’s tidal radius is that the ma ss ratio of the secondary t o 
the primary be less than ~0.33 llWhitehurst fc Kindi 19911 ) . 
If we assume that the secondary is a main sequence star that 
fills its Roche lobe, then the relation m 2 — 0.11 Phr holds, 
where m 2 is the mass of the secondary in solar u n its a nd 
Ph r is the orbital period in hours fe.g. lFrank. King fc Raiiiel 
EDift. Hence any accreting binary that has a short enough 
orbital period may exhibit superhumps. Indeed, there exists 
a class of short-period, p ersistently bright CYs that ex hibit 
perm anent superhumps ( Pattersonl Il99ij ; llRetter fc NavloJI 

IH as well et alJ (1200111 discuss analogous superhump be¬ 
haviour in LMXBs. Although superhumps have been found 
in the optical lightcurves of several black hole and neu¬ 
tron star LMXBs, they cannot be produced by the same 
mechanism, since the optical output of X-ray brig ht LMXB s 
is dominated by reprocessed X-rays. Instead. lHaswell et all 
( 2001 ) proposed that the modulation is due to variation of 
the solid angle that the disc subtends to the X-ray source 
(on the superhump period); in this model, superhumping 
might be expected to be more prominent in low inclination 
systems. 

4U1916—053 is a neutron star LMXB, with an X- 


Observation 

Date 

Exp 

Filter 

Al 

2000 Jun 25 

34 ks 

Medium 

A2 

2001 Jun 29 

56 ks 

Medium 

A3 

2002 Jan 6 

61 ks 

Thin 

P4 

2004 Jul 17 

18 ks 

Medium 

P5 

2004 Jul 19a 

22 ks 

Medium 

P6 

2004 Jul 19b 

27 ks 

Medium 


ray period of 50.00 min and an optical period of 50.4589 
min dCallanan et alJll995l) . lHaswell et all ll200ll) show that 
LMXBs with orbital periods shorter than ~4.2 hr are likely 
to exhibit superhumps, and identified 4U1916—053 as a 
persistent irradiated superhumping source. 4U 1916—053 is 
a high inclination s ystem, with periodic i ntensity dips in 
the X-ray lightcurve dWhite fc SwanlJll982i) . These dips are 
due to photo-electric absorption by material on the outer 
edge of the accretion disc; many believe that the inflated 
bulge on the outer disc rim that is caused by the colli¬ 
sion betwee n the gas stream and the outer disc is responsi¬ 
ble (see e.g. IWhite fc IIcdTlri f)82l) . The X-ray modula tion of 
4U 1916—053 shows striking variability dSmale et alJIlOSsl) . 
These variations repeat on a ~4 day period, an d are caused 
by th e prece ssion of the a c cretio n disc (see e.g. IChou et alJ 
1200ll) . When lRetter et ‘ahl d2002t) made power density spec¬ 
tra (PDS) of the X-ray lightcurve of 4U 1916-05, they found 
peaks corresponding to both the X-ray and optical periods. 
After removing the dipping intervals, the optical peak was 
removed; hence the di ps we re shown to occur on the super¬ 
hump period. iRetter et all (’ 20021) concluded that that the 
observed superhumps arose from the same thickened region 
of the outer disc that caused the absorption of the X-rays, 
allowing superhumps to be seen in high inclination LMXBs. 


3 OBSERVATIONS AND DATA ANALYSIS 

In additi on to the three XM M-Newton observations anal¬ 
ysed by iTrudolvubov et alJ ll2002fi . we conducted a pro¬ 
gramme of four ~20 ks observations over 2004, July 16-19. 
We present results from our analysis of the archival data, 
along with three of the four 2004 observations; the other 
observation suffered from flaring in the particle background 
over 90% of the observation and is not considered further 
here. A journal of observations is presented in Table Q 

We analysed data from the pn, MOS1 and MOS2 in¬ 
struments, which share the same 30” x 30” field of view. We 
used version 6.0.0 of the SAS software suite 1 to obtain the 
data products, as well as the latest calibration data. For each 
observation, we selected a circular extraction region with a 
40” radius, centred on Bo 158, and an equivalent source-free 
region for the background. The background region was on 


1 http://xmm.vilspa.esa.es 
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Figure 1. Combined EPIC 0.3—10 keV lightcurves (in count s — *) 
of Bo 158 from the archival XMM-Newton observations, A1-A3; 
the dates of each observation are shown. 

the same chip as the source, and at a similar angular off¬ 
set from the optical axis. We extracted lightcurves from the 
source and background regions in the 0.3-10 keV, 0.3-2.5 
keV and 2.5-10 keV energy bands with 2.6 s binning; these 
were analysed using FTOOLS version 5.3.1. We also ob¬ 
tained pn spectra of the source and background regions in 
4096 channels of 5 eV, and generated response matrix and 
ancillary response files from the source spectral files. The 
spectra were then grouped for a minimum of 50 counts per 
bin. Spectral analysis was performed using X'SPEC 11.3.1. 


4 RESULTS 

The 0.3-10 keV EPIC (MOS1 + MOS2 + pn) lightcurves of 
Observations A1-A3 are presented in Fig.0 with x- and y- 
axes set to the same scale, and with 200 s binning. Most 
striking is Observation A3, with six dipping interva ls on 
a 10017T50 second period (ITrudolvubov et al.l l2003) : the 
structure of the dipping is seen to vary substantially even 
over one observatio n of 60 ks. We emphasise that the 30% 
dipping reported bv lTrudolvubov et alJ (l2002fi for observa¬ 
tion A1 is an overall average; a deep dip is seen at ~28 ks 
into the observation, but very little evidence of dipping is 
observed in the intervals of expected dipping at ~8 ks and 
~18 ks into the observation. In this regard, Observation A1 
resembles the 1985, October lightcurve of 4U 1916—053, ob¬ 
served by EXOSAT, where deep dips were observ ed only 
after the first four orbital cycles dSmale et a nma. Little 
evidence of variability is seen in the lightcurve of Observa¬ 
tion A2. 

In Fig. |21 we present the 0.3-10 keV EPIC lightcurves 
of Observations P4-P6; the x-axis is scaled to P6, and the 
y-axis matches that of Fig.0 The most prominent feature is 
the dip in P5; it has a depth of ~100% and a total duratio n 
of ~2500 s. Using the period of iTrudolvubov et alJ ll2002lb 
we identified the expected times of dipping, labelled ‘D’, in 
P4 and P6, using the deepest part of the dip in P5 as phase 
zero. 

In the P4 lightcurve, there is no evidence for dipping 
during the first expected dip interval, but some evidence of 



Figure 2. Combined EPIC 0.3—10 keV lightcurves (in count s —1 ) 
of Bo 158 from the 2004 XMM-Newton observations, P4-P6; the 
dates of each observation are shown. The horizontal line in the 
P5 lightcurve indicates the interval used for spectral fitting. The 
shaded area of the P6 lightcurve indicates a period of background 
flaring. Times when dipping is expected are labelled “D”, using 
the deepest part, of thedip in P5 as time zero, and the period of 
|TUidolYubov_et_alJ j2002|). 

dipping ~4000 s after the second interval, 20 cycles away 
from the dip in P5. Hence, this possible dip would require a 
period that is eit her ~200 s shorter or ^320 s longer than the 
10017 s given bv lTrudolvnhov et al ] (l2002i) . However, there 
is no other evidence for these other periods in the lightcurves 
of P4, P5 or P6; hence it is clear that the dipping behaviour 
of Bo 158 evolves on a time scale of a few days, just like that 
of 4U 1916-053. 

Several emission models were fitted to the 0.3-10 keV 
pn spectrum of P4, each suffering absorption by material 
in the line of sight. P4 was chosen because it had the 
longest interval of persistent emission that was not contam¬ 
inated by background flares; the resulting source spectrum 
contained ~2200 counts. We applied the two models that 
ITrudolvubov et al.l (I2002T) used to model the spectra of Al— 
A3, namely a power law model and a Comptonisation model 
(comptt in XSPEC). We also applied a two component 
model. The emission of many Galactic LMXBs has been 
succesfully described by a model consisti ng of a black body 
and a cut-off power law fe.g. Ichurch fe Baluci nska-Churcbl 
Il995l : IChurch et al.lll993 iBarnar^^l^d'll^OO^i : we approxi- 
mate this model to a blackbody + power law model, because 
of the narrow pass-band. Table |2]shows the best fits to the 
spectrum with each model; uncertainties are quoted at the 
90% confidence level. 

We find that the best fit parameters for the power law 
an d COMPTT models a gree w ell with the values presented 
bv iTrudolvubov et alJ ( 2002 1: however, the two-component 
model provided the best fit. We present the unfolded SED 
for the two component fit in Fig. [3 the emission from the 
blackbody (BB) and power law (PO) components are shown 
separately. We see that the blackbody dominates the emis¬ 
sion above 1.5 keV. We find a 0.3-10 keV flux of ~2xl0 -12 
erg cm -2 s -1 for all the fits to the P4 data; this gives a 
0.3-10 keV luminos ity of ~1.4x l0 38 erg s -1 , assuming a 
distance of 760 kpc llvan den BerglfeoPOTl . 
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Table 2. Best fits to the P4 (2004 July 16) pn spectrum of Bo 158 for several emission models; these are power law (po), Comptonisation 
(COMPTt), and blackbody + power law (bb+po) models, suffering line-of-sight absorption by the interstellar medium. Ah is the equivalent 
hydrogen column density of the absorber, To is the seed photon temperature (COMPTT), T is the blackbody temperature (COMPTT, bb+po), 
r is the optical depth (comptt) and Y is the photon index (po, bb+po). The goodness of fit is given by x 2 /dof. and To. 3—10 is the 
unabsorbed flux in the 0.3—10 keV band. 


Model 

A h 

atom cm -2 

kT 0 

keV 

kT 

keV 

T 

r 

X 2 /dof 

To. 3-10 

10 -12 erg cm 2 s -1 

PO 

0.1 




0.57+0.09 

30/24 

2.3+0.3 

COMPTT 

0.1 

0.08+0.07 

2.4+0.6 

17+2 


22/19 

o 9+0.2 

z.z_2 0 

BB+PO 

0.1 


2.0+0.2 


2.0+0.3 

19/19 

2.1+0.8 


We then extracted the SED from the interval of 
^persistent emission in P5 indicated in Fig. 3 by a horizontal 
line. The resulting background-subtracted SED contained 
452 counts in the 0.3-10 keV band, which we divided into 
10 spectral bins. The spectral shape of the P5 SED was con¬ 
sistent with that of P4, with a luminosity of 1.2±0.2xl0 38 
erg s- 1 . Hence the luminosity of persistent emission in P5 
is consistent with that of P4 within 3 a. 

The depth of dipping in Al varies from ~0 to ^70% 
with no significant change in the mean intensity, suggesting 
that the amplitude of dipping is not simply anticorrelated 
with the source luminosity. Instead, the variation in dipping 
behaviour may be caused by disc precession. This hypothesis 
motivated our simulation of the accretion disc in Bo 158, 
using three dimensional smoothed partical hydrodynamics, 
discussed in Sect.[o] 

The lightcurves of XMM-Newton observations of Bo 158 
show no eclipses; hence, we know that we are not viewing 
the system edge on. If the disc were tilted with respect to 
the binary plane, and precessing, then one might expect to 
observe dips in some part of the disc precession cycle, but not 
in others. Dipping is observed throughout observation A3; 
this suggests that the dipping phase in the disc precession 
cycle lasts >60 ks. The Al lightcurve covered three intervals 
of expected dipping, yet only one dip is seen, toward the 
end of the observation; we suggest that this dip signals the 
onset of the dipping phase. Contrariwise, P5 and P6 appear 
to sample the end of the dipping phase, as a dip is observed 
in P5, yet no dips are seen in P6. 

5 SPH SIMULATION OF THE DISC 
5.1 Binary Parameters 

The accretion disc was modelled using a three-dimensional 
Smoothed Particle Hydrodynamics (SPH) co m put e r cod e 
that has been des cribe d in detail in [Murrayl fl99f , H198), 
iTruss et alJ ll200fl) , and iFoulkes. Haswel^^Turrav l 200, c ). 
We assumed the orbital peri od to be 10017 s, as obtained 
bv iTrudolvubov et, alJ d2002l) . Table 3 lists the system pa¬ 
rameters adopted. 

The dipping source Bo 158 is a bright globular cluster 
X-ray source, with a 2.78 hr binary period. Thirteen Galac¬ 
tic globular clusters contain bright X-ray sources; twelve of 
these are neutron star LMX Bs, while the primary of the 
other one is unknown (see e.g. lin’t Zand et all2004i) . Hence 


Bo 158 is a likely neutron star LMXB, and we assume the 
primary mass to be 1.4 Mg. 

For the secondary, we considered a main sequence star 
and a white dwarf, s ince 4U 191 6—053 has a likely white 
dwarf secondary (e.g. IChou et al ] 2001). For a Roche lobe- 
filling main seq uence star t h e appr oximate relation m ~ 0.11 
Phr holds (e.g. iFrank et alJ 1200 2). giving a mass of ~0.30 
Mg. If instead the star is a white dwarf, using the mass 
radius relation of iNauenberd (1972) and Roche geometry 
gives an implausibly small mass of 0.005 Mg. Assuming a 
main sequence secondary, we found the mass ratio to be 0.2, 
indicating that superhumps and disc precession were likely. 

Finally, the luminosity of the system was taken to be 
1.4xl0 38 erg s _1 , the 0.3-10 keV luminosity of Bo 158 in 
Observation P4. Such a high luminosity may be expected 
to cause warpin g of the accret ion disc, even for a previously 
flat disc (see e.g . lPringleillf996 ); warping is discussed in Sect. 
5.3. 

The accretion disc had an open inner boundary con¬ 
dition in the form of a hole of radius ri = 0.025a, where 
a is the binary separation, centred on the position of the 
primary object. Particles entering the hole were removed 
from the simulation. Particles that re-entered the secondary 
Roche lobe were also removed from the simulation as were 
particles that were ejected from the disc at a distance > 0.9a 
from the centre of mass. 

We assumed an isothermal equation of state and that 
the dissipated energy was radiated from the point at 
which it was genera t ed, as electromagnetic radiation. The 
IShakura fc Sunvaevl (1973) viscosity parameters were set 
to aiow = 0.1 and ahigh = 1.0, and the vi scosity state 
changed smoothly as described in ITruss et alJ 1200(1) . The 
SPH smoothing length, h, was allowed to vary in both space 
and time and had a maximum value of 0.01a. 


5.1.1 The gas stream 

We simulated the mass loss from the secondary by introduc¬ 
ing particles at the inner Lagrangian point (Ti). The mass 
transfer rate and the particle transfer rate were provided as 
input parameters, and the mass of each particle was derived 
from these parameters. A particle was inserted with an ini¬ 
tial velocity (in the orbital plane) equal to the local sound 
speed of the donor, cd, in a direction prograde of the binary 
axis. The z velocity of the inserted particle was chosen from 
a Gaussian distribution, with a zero mean and a variance of 
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unfolded spectrum 



Energy (keV) 


Figure 3. Unfolded EPIC-pn SED of Bo 158 from observation 
P4, described by the best fit two component model. The dotted 
line represents the blackbody component (BB), while the dashed 
line shows the power law component (PO). The solid line shows 
the sum of the two components. The blackbody component dom¬ 
inates the flux above 1.5 keV. 

O.Icd. However, the inflation of the site of collision between 
the gas stream and outer disc was not modelled. 

5.1.2 The initial non-warped accretion disc 


rithm was used to determine regions of self-shadow. For each 
particle found on the disc surface a light-ray was projected 
from the particle to the position of the radiation source at 
the centre of the disc. The particle was deemed to be il¬ 
luminated by the radiation source if this light-ray did not 
intersect any disc material between the particle surface po¬ 
sition and the radiation source (i.e. the particle could see 
the central radiation source). The radiation force was only 
applied to particles that were considered to form part of the 
disc surface and were illuminated by the central radiation 
source. 

5.3 Disc warping and precession measure 

For an optically thick disc, a warp can develop as a result 
of the radiation force (e .g. IPrin [Qgilvie fc Dubuj 

120011 iFoulkes et aljEooijl) . This is due to the fact that any 
radiation absorbed by a specific region on the disc surface 
will be later re-radiated from the same spot, normal to the 
disc surface. Hence any anisotropy in the disc structure will 
cause an uneven distribution of back-reaction forces on the 
disc surface, further perturbing the disc. A sufficiently high 
luminosity can in duce and s ustain a warp even in an origi¬ 
nally flat disc llPrin gle 1996, and r eferences t herei n). 

_The two measures defined by lLarwood fc Papaloizoul 

(1997) were used to measure the disc warping and the 
amount of warp precession. They defined an angle j as the 
angle between the total disc angular momentum vector and 
the angular momentum vector for a specific disc annulus, 
i.e. 


The simulation was started with zero mass in the accretion 
disc and with the central radiation source switched off. A sin¬ 
gle particle was injected into the simulation every O.Olfl”^ 
at the L\ point as described above until a quasi-steady mass 
equilibrium was reached within the disc. This was taken to 
be when the number of particles inserted at the L\ point, 
the mass transfer rate, was approximately equal to number 
of particles leaving the simulation at the accretor, the ac¬ 
cretion rate. The simulations were continued for another 3 
orbital periods to ensure mass equilibrium. The number of 
particles in the simulated accretion disc was approximately 
40,000 giving a good spatial resolution; the average number 
of ‘neighbours’, i.e. the average number of particles used in 
the SPH update equations, was 8.2 particles. The simulated 
disc encountered the Lindblad 3:1 resonance and became ec¬ 
centric. The disc precessed in a prograde direction giving rise 
to superh umps in the simulated dissipation light-curves (c.f. 
IFoulkes et al.lli2004 '). The radiation source was then turned 
on which gave rise to a very small number of particles being 
ejected from the accretion disc. 


5.2 Surface finding algorithm & self-shadowing 

Accretion-powered radiation from the inner regions of the 
disc and the accreting object itself exer t a for ce on the irra¬ 
diated disc surface. Following IPrin iH ( 19961 1 the radiation 
source is modelled as a point source at the centre of mass 
of the accretor. To apply this force, particles on the surface 
of the accretion disc had to be identified. We used a convex 
hull algorithm to find the surface part icles as described in 
iMurravl il998fl and lFoulkes et al . ( 12005 1. A ray-tracing algo¬ 


J A ■ 3d 

Jj4 | | 3 d 

The term 3a is the total angular momentum within the 
specific annulus and was calculated by summing the angular 
momentum for each particle within the annulus. The term 
Jd is the total disc angular momentum and was calculated 
by summing all the angular momenta for all particles within 
the disc. An angle n was also defined which measures the 
amount of precession of the disc angular momentum relative 
to the initial binary orbital angular momentum, Jo 


cos n = 


(Jo x 3d) ■ u 
I Jo X Jo 11 u I 


( 2 ) 


where u is any arbitrary vector in the binary orbital plane. 


5.4 Numerical modelling results 

Prior to irradiation, the disc was asymmetric about the bi¬ 
nary axis and precessed in a prograde direction relative to 
the inertial frame. After switching the radiation source on 
we ran the model for a further 50 orbital periods, and found 
that this illumination introduced a warp in the disc. When 
the radiation source was removed, the warp would dissipate 
and the disc would return to the orbital plane. 

As a result of the disc precession, viscous stresses in the 
disc vary significantly with time. Figure 0 shows how the 
resultant energy dissipation in different regions of the disc 
varies with time. The disc luminosity was not modelled in 
detail. We assume that the luminosity was directly related 
to the disc regions with significant energy release through 
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Table 3. Binary system parameters for the system modelled. The columns are: the total system mass, the system orbital period, the 
rate of mass loss from the secondary, the system mass ratio, the physical luminosity of the central radiation source and the ratio of the 
physical luminosity to the Eddington limit for the system. 


Parameter Mt 

Mq 

Porb 

days 

Msec 

M^yr- 1 

q 

(Ma/Mi) 

Physical luminosity ( L *) 
(erg s _1 ) 

Z/* / L e dd 

Value 1.8 

0.1159 

2.22 X 10“ 8 

0.2 

1.4 X 10 38 

0.8 


viscous dissipation. The viscous dissipation heats the gas in 
the accretion disc and it is assumed that the heat is radiated 
away from the point at which it was generated. Superposed 
on a steady signal there is a repeating series of “humps.” 
The spacing of the humps corresponds to the superhump 
period. The humps consist of three separate major compo¬ 
nents. The maximum occurs at phase ~0.2, followed by a 
secondary maximum at phase ~0.5, and the minimum at 
^0.7; however the relative strengths of the humps vary from 
cycle to cycle. We stress that the dissipation lightcurves are 
only a diagnostic for the disc properties, i.e. for identifying 
the the superhumps, as well as the intervals when the disc 
structure is most extended or most compressed. They are 
not representative of true optical lightcurves. 

In order to determine the superhump period, P s h, we 
obtained a power density spectrum from ~30 superhump 
cycles of the simulated light curve. We estimated the su¬ 
perhump period to be (1.035 ± 0.005)P or b- This implies the 
precession period of the outer regions, _P prec = (29 ± l)P 0 rb, 
or 81±3 hr. 

Fig . [^contains projection plots for the time period cor¬ 
responding to the peak labelled “maxi” in Fig. Q] at phase 
~2.2. A very strong spiral density compression wave can 
be seen at the upper edge of the disc. This wave is so in¬ 
tense that it is removing material from the accretion disc 
an d returning it_ back to the Roche lobe of the secondary, 
see iFoulkes et al* 1 2004 1 for a full detailed description of a 
similar system with a mass ratio of 0.1. 

The two upper right-hand plots of Fig. 0 labelled yz- 
view and xz-view, are side views of the disc in the y-z and x-z 
directions respectively. The yz-view plot is a projection view 
of the disc as seen from the secondary, similarly the xz-view 
is a projection plot with the secondary located to the right 
of the plot. The disc warp is clearly apparent in these two 
plots. The warp is odd symmetrical about the centre of the 
disc. The maximum value of the warp is located at a distance 
approximately 0.1a either side of the primary position, see 
yz-view of Fig.^1 

The lower plot of Fig. [£] labelled accretor-view, shows 
the distribution of the particles as seen from the compact 
object. The horizontal axis is the binary orbital phase, L i is 
located at phase 0 and disc material flows from right to left 
with the stream-disc impact region located at approximately 
phase 0.9. The vertical axis is the elevation angle of the 
particle as seen from the primary position. From this plot it 
can be seen that the radiation force has pushed disc material 
out of the orbital plane. The warp reaches a maximum height 
above the orbital plane at phases ~0.1 and ~0.6, and has 
minima at phases ~0.45 and ^0.95. We also see that the 
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Figure 4. The top left plot shows SPH viscous dissipation light 
curves for different regions of the accretion disc. The top curve, 
labelled a, is for the whole accretion disc, then each curve in 
descending order corresponds to dissipation in the disc at radii 
> 0.05a(b), > 0.1a(c), > 0.2 a(d) and > 0.3a(e) respectively. The 
light curve minimum and two maximum points are indicated. The 
lower left-hand plot is the signal from the disc inner region and 
was generated by subtracting light curve b from light curve a. 
The plots on the right-hand side are disc dissipation maps that 
correspond to the light curve minimum and two maxima. The 
circles labelled b, c, d and e show the distances from the primarily 
that correspond to the light curves b, c, d, e respectively. 


disc remains mainly in the orbital plane, although it can be 
seen that there is a small S-wave in the structure of the disc. 

The warp amplitude and size precessed as a solid body 
in a retrograde direction relative to the inertial frame. The 
warp precession period, Pwarp, was determined by applying 
equation © to each time step of the simulation. This gave 
a precession angle relative to some arbitrary start angle for 
each simulation time step. The warp precession rate was 
then determined by fitting a straight line to th ese data us- 
ing a Numerical Recipes least squares method, llPress et alJ 
I198 (tTi . The precession rate was found using the gradient of 
the line extracted from the least squares fit. We found that 
-Pwarp ~11 Porb- FigureH3shows the radial profile of the warp 
for five consecutive orbital cycles; the maximum extent of 
the warp is 6-11° above the plane of the disc. 
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Figure 5. Particle projection plots for the SPH model. The position of each particle is indicated by a small black dot. The plot labelled 
xy-view is a plan view of the accretion disc as seen from above the disc. The cross at the centre of the plot shows the position of the 
primary object. The solid dark line is the Roche lobe of the primary and the Li point is to the right and middle of the plot. The two 
plots xz-view and yz-view are particle projection plots on a plane perpendicular to the orbital plane and through the system axis. The 
bottom plot, accretor-view, shows the particle distribution as seen from the compact object. The horizontal axis is the orbital phase, 
the Li point is at phase 0 and the stream/disc impact region is at approximately phase 0.9. The vertical axis is the angle, in degrees, 
between a particle and the orbital plane when viewed from the compact object. The disc material flows from right to left. 


6 DISCUSSION 

We report a single ~100% dip in the 0.3-10 keV EPIC 
lightcurve of Bo 158 from Observation P5, but find lit- 
tle evidence of d i pping in Observations P4 and P6. 
iTrndolviihov et~, §2QQ2) propose two models to explain 
the energy independence of dipping that they inferred from 
Bo 158: (1) absorption by a highly ionised region and (2) 
partial covering of an extended source by an opaque ab¬ 
sorber that occults varying fractions of the source. However, 
in their figure showing both non-dip and dip SEDs, the two 
SEDs converge at high energies, showing that the dipping is 
indeed energy dependent. 

The dipping behaviour of many Galactic LMXBs has 
been well described by a single model: absorption of a point¬ 
like blackbody plus progressive covering of an extend ed 
emission region by an extended absorber llChurch et alJ 


Il997t) . This extended emission is caused by unsaturated 
inverse-Comptonisation of cool photons on hot electrons 
in an accret ion disc corona that has a radius o f ~10,000- 
500,000 km ( IChurch fc Balucihska-Churchlll2004ll . The satu¬ 
rated dipping exhibited by 4U 1624—490 during a 1985 EX- 
OSAT observatio n is particula rly strong evi den ce for an ex¬ 
tended absorber llChurch fc Bahiciriska-Churcblfl995 '). The 
0.1-200 keV observations of 4U 1916—053 with Beppo-SAX 
have shown that the absorber is so dense that ~100% photo¬ 
electric absorption o ccurs up to 10 keV; in fact, the dipping 
is seen up to 40 keV llChurch et all 199811 . However, the high 
luminosity of Bo 158, together with the large blackbody con¬ 
tribution, mean that the dipping is unlikely to be energy- 
independent. 

Disc precession is inferred from the 0.3-10 keV 
lightcurves of Bo 158, as is expected given its extreme mass 
ratio (short orbital period). As such, it resembles the Galac- 
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Figure 6. Radial warp profiles evaluated using equation 0 . 
The vertical axis is the warp amplitude. The horizontal axis is 
distance from the primary object normalized such that the bi¬ 
nary separation is 1. Plots (a), (b), (c), (d) and (e) are for five 
consecutive orbital cycles. 

tic superhumping LMXB 4U 1916—053. Since the LMXB Bo 
158 is in a globular cluster near the centre of M31, it is un¬ 
likely that the optical period will ever be known. However, 
our Fourier analysis of the simulated dissipation lightcurves 
indicates a superhump period that is 3.5±0.5% longer than 
the orbital period. Given the asso ciation betwe en the dips 
and superhump period reported bv lRetter et al. ( 12002 1. the 
10017 s period may be the superhump period, in which case, 
the orbital period would be ~4% shorter. Such shortening 
of the period would not dramatically affect the outcome of 
our SPH modelling. 

Our simulations of the disc show two distinct types of 
variability in the disc structure. First is the elongation and 
prograde precession of the disc due to tidal interactions with 
the secondary at the 3:1 resonance; the disc precesses on 
period of 81±3 hr. We also see warping of the accretion disc, 
driven by irradiation of the disc surface by the central X-ray 
source; the warp is stable exhibits retrograde precession on 
a ~31-hr period. 

It is therefore important to establish which region is re¬ 
sponsible for the observed variation in dip morphology. The 
lightcurves of observations A1-P6 show no eclipses. From 
Kepler’s law and th e rat io of the secondary ratio to the bi¬ 
nary separation ( Eggletonlll983il . the secondary has an an¬ 
gular radius of ~15°; hence the inclination <75°. We see 
from Fig. 6 that the disc warp does not deviate from the 
plane of the disc by more than 11° in our simulations, sug¬ 
gesting that the observed dips are likely to evolve on the 
disc precession period. 


7 CONCLUSIONS 

We have analysed three new XMM-Newton observations of 
the M31 dipping LMXB Bo 158, in addition t o re-a nalysing 
the t hree observations discussed in iTrudolvubov et all 
i2002h . The newer observations spanned ~3 days in 2004, 
July. We find that that the relationship between source in¬ 


tensity and depth of dipp ing is not simple, as described by 
ITrudolvubov et~ ( 12002 1. Instead, we believe that the ob¬ 
served variation in dipping behaviour is caused by precession 
in the accretion disc; dipping would be confined to a limited 
phase range in the disc precession cycle. 

We modelled the accretion disc with 3D SPH, and found 
prograde disc precession on a 81±3 hr period, as well as 
radiatively driven disc warp that precessed on a 31 hr period 
in a retrograde fashion. We find that the disc precession is 
most likely to affect the observed dipping behaviour. Hence, 
we predict that the dipping behaviour of Bo 158 experiences 
a 81±3 hour cycle; this period is consistent with the observed 
variation of the dips. 
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